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1.  INTRODUCTION 


This  report  documents  the  use  of  TMR  processing  techniques  to  produce 
GSRS  spin  and  coning  motion  measurements  from  data  collected  by  the  instru¬ 
mentation  radars  at  WSMR.  These  measurements  are  made  from  single-site 
data.  Two  software  packages  have  been  developed  to  make  the  measurements. 
Program  WSGSRS  measures  spin  frequency,  and  program  TCM  measures  coning 
motion.  These  programs  have  been  implemented  on  both  the  Univac  1108  and 

PDP-11  Radar  Graphics  Laboratory  at  WSMR. 

Some  familiarity  with  TMR  processing  on  the  part  of  the  reader  is  as¬ 
sumed.  This  familiarity  could  be  gained  by  a  perusal  of  the  following 
reports. 

Radar  Measurement  Performance  Improvement  by  Motion  Resolution 
Methods  (U),  A.  W.  Rihaczek,  W.  B.  Kendall,  G.  W.  Lank,  MARK 
Resources  Report  110-11,  Final  Report  under  Contract  DAAD07- 

75- C-0046,  June  1975,  AD-C002  333,  CONFIDENTIAL. 

Development  of  Improved  Data  Processing  Methods  (U) ,  A.  W.  Rihaczek, 

P.  A.  Gartenberg,  W.  B.  Kendall,  J.  J.  Heimhold,  I.  P.  Bottlik, 

MARK  Resrouces  Report  123-14,  Final  Report  under  Contract  DAAD07- 

76- C-0019,  May  1976,  AD-C006  470,  CONFIDENTIAL. 

Feasibility  Investigation  of  Motion-Resolution  Processing  for 
Artillery  Shells  (U) ,  A.  W.  Rihaczek,  W.  B.  Kendall,  J.  L.  Dessinger, 
MARK  Resources  Report  123-23,  Final  Report  under  Contract  DAAD07- 
76-C-0019,  Modification  P00001,  December  1976,  AD-C009  233, 
CONFIDENTIAL. 

Improved  Processing  for  Instrument at ion- Radar  Data  (U) ,  W.  B.  Kendall, 
A.  W.  Rihaczek,  P.  A.  Gartenberg,  G.  W.  Lank,  J.  L.  Dessinger,  MARK 
Resources  Report  137-36,  Final  Report  under  Contract  DAAD07-77-C- 
0028,  December  1977,  AD-C013  495,  CONFIDENTIAL. 

Familiarity  with  Section  3.2  of  the  last  listed  report  is  especially  im¬ 
portant  . 
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2 .  DATA  REDUCTION 
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Before  eargec  morion  measurements  can  be  made,  the  radar  data  on  Che 
WSMR  FPS-16  or  MPS-36  site  tape  must  be  reduced  via  the  WSPRC  software. 

Both  the  spin  frequency  and  coning  motion  programs  use  data  which  have 
been  compensated  for  the  gross  trajectory  motion  of  the  missile.  The  pro¬ 
cessing  required  by  Che  spin-frequency  program  will  be  described  first. 

In  order  to  measure  spin  frequency,  the  data  must  be  compensated  to 
clearly  show  spin  lines  on  a  TMR  plot  such  as  Figure  1.  The  spin  lines 
will  give  reliable  measurements  only  when  Che  central  peak  representing 
the  base  return  is  well  aligned  at  zero  Doppler.  The  compensation  (for 
the  trajectory)  which  aligns  Che  central  peak  is  carried  out  in  two  steps: 
(1)  Che  initial  correction;  and  (2)  the  lowpass-f ilter  correction. 

2 . 1  INITIAL  CORRECTION 

The  result  of  the  initial  correction  is  compensated  data  for  which  the 
central  spectral  peak  (from  the  missile  base)  remains  between  +  20  Hz.  To 
produce  Che  initial  correction,  Che  site  tape  is  processed  by  the  WSPRC 
software  with  the  following  parameters  specified: 

SPARMS 

LEVEL-1,1,1,1,1, 

GENT AP-. TRUE. , 

TB IAS- (launch  t ime ) , 

PRF-320.1729, 

FREQ- (frequency) , 

AUT0R-.TRUE. , 

XKNTSR-0. , 

SUBDVR-4, 

NAVGR-8 , 

SUBDVP-3 , 

NAVGP-8 , 

BRKMNP-0. , 

D0PPLT-1 . , 

NFFT-256, 

LENWDW»64 , 

LAG-16, 

SEND 

Figure  2  is  an  example  of  a  satisfactory  initial  correction  produced 
by  such  a  run  of  WSPRC.  When  the  initial  fit  satisfactorily  aligns  the 
central  peak  as  in  Figure  2,  processing  continues  with  the  second  step, 
the  lowpass-f ilter  correction. 
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Figure  1.  TMR.  Plot  Showing  Spin  Lines  at  Multiples  of 
the  Spin  Frequency 
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Occasionally  the  initial  run  of  WSPRC  will  not  produce  satisfactory  re 
suits.  The  TMR  plots  may  show  sections  of  bad  data  or  clutter.  Figure  3 
shows  the  result  of  an  initial  correction  applied  to  data  for  which  the  re¬ 
corded  range  was  invalid  between  9.5  and  11.5  seconds.  In  the  bad  section 
of  data  the  central  peak  wanders  all  over  the  page.  In  such  a  case  it  is 
best  to  rerun  the  WSPRC  initial  correction,  setting  parameter  EDIT=5,  and 
using  edit  cards  to  specify  the  time  interval  for  which  the  range  data  are 
bad.  Figure  4  shows  the  result  of  the  WSPRC  initial  correction  for  which 
bad-data  interval  is  specified  as  9.5  to  11.5  seconds.  The  main  peak  now 
oscillates  about  zero  Doppler,  but  stays  within  +  50  Hz . 

When  the  initial  correction  does  not  center  the  central  peak  between 
+  20  Hz,  a  spline  may  be  generated  by  hand  to  improve  the  correction.  The 
R-dot  Peaks  (or  Doppler  Peaks)  plot  generated  by  setting  parameter  D0PPLT=1 
in  LEVEL (5)  of  WSPRC  is  useful  in  constructing  the  hand  spline.  Figure  5 
shows  the  R-dot  Peaks  plot  for  the  data  of  Figure  4.  The  first  step  in 
creating  the  hand  spline  is  to  select  the  breakpoints  at  the  peaks  and 
valleys  of  the  R-dot  Peaks  plot,  as  shown  in  Figure  6.  The  first  few  break 
points  are: 


BREAKPOINT  TIME 

1  .  4.82 

2  5.61 

3  6.50 

4  6.92 

5  7.35 

Next,  the  values  of  R  at  the  breakpoints  are  read  from  the  plot: 


BREAKPOINT 

TIME 

R  (m/sec) 

1 

4.82 

.05 

2 

5.61 

-.02 

3 

6.50 

.02 

4 

6.92 

-.06 

5 

7.35 

.09 

10 


Figure  6.  Picking  the  Breakpoints  and  Drawing  the  Spline 
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Finally,  the  range  is  computed  at  each  breakpoint,  and  the  spline  is  com¬ 
plete.  The  range  value  at  the  right  of  each  breakpoint  internal  is  given 
by  a  formula  which  uses  the  previously  computed  range  value  at  the  left 
of  the  interval,  and  the  time  and  R  values  at  the  left  and  right  of  the  in¬ 
terval.  The  formula  is 

\  ■  \  +  7  (l W  (VrL>  • 

Zero  is  used  for  the  first  range  value,  and  the  succeeding  range  values 
are  calculated  according  to  the  recursive  formula.  In  this  example  the 
second  range  value  is  given  by 

Rn  -  0.  +  |  ( . 05+(-. 02) ) (5.61-4.82)  =  .01185  m  . 

This  value  becomes  for  the  computation  of  the  next  range  value: 

Rr  *  .01185  +  j  (-.02+. 02) (6.50-5.61)  -  .01185  m  . 

Continuing  this  process,  the  first  five  points  of  the  spline  are: 


BREAKPOINT 

TIME 

R(m) 

R  (m/sec) 

1 

4.82 

0. 

.05 

2 

5.61 

.  011*85 

-.02 

3 

6.50 

.01185 

.02 

4 

6.92 

.00345 

-.06 

5 

7.35 

.0099 

.09 

These  spline  data  are  entered  into  the  computer,  and  then  read  in  LEVEL (4)  of 
WSPRC  by  setting  parameters  104*5  (for  cards)  and  105=0  in  namelist  $PARMS. 
Figure  7  shows  the  result  of  applying  the  spline  correction  from  Figure  6  to 
the  data  in  Figure  4.  One  more  hand  spline  correction  is  necessary  to  further 
correct  the  data  before  proceeding  to  the  lowpass  filter  correction. 
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2.1.2  Processing  with  Ground  Clutter 


Another  problem  which  may  appear  during  the  initial  correction  is  ground 
clutter.  Figure  8  shows  ground  clutter,  the  strong  line  moving  from  the 
bottom  of  the  plot  to  the  upper  right-hand  corner.  As  long  as  the  ground 
clutter  does  not  become  confused  with  the  central  peak,  correction  of  the 
data  may  continue  as  usual.  The  ground  clutter  will  become  a  nuisance 
later,  though,  because  it  will  interfere  with  the  spin  lines. 

2.1.3  Missing  Data 

Occasionally,  an  interval  of  data  may  be  no  good  at  all,  and  impossible 
to  correct.  Such  a  case  is  shown  in  Figure  9  between  the  times  24  and  26 
seconds.  In  that  figure  the  central  peak  from  the  missile  base  seems  to  dis¬ 
appear.  Lowpass  filtering  fails  to  provide  any  clues  as  to  what  was  happening 
here,  so  this  section  of  the  data  should  be  ignored.  Most  often,  it  will  be 
possible  to  correct  the  data  fairly  well  by  using  splines  wherever  necessary. 

2.2  LOWPASS  FILTER  CORRECTION 

Once  the  initial  correction  is  good  enough  so  that  the  central  peak 
stays  between  -20  and  +  20  Hz,  the  data  may  be  further  corrected  by  lowpass 
filtering  that  and  using  the  lowpass  filter  (LPF)  data  as  input  to  LEVEL (3) 
of  WSPRC. 

LEVEL (3)  fits  a  spline  to  the  LPF  data,  and  the  spline  is  used  in 
LEVEL (4)  to  further  correct  the  initially-corrected  data.  To  do  this,  the 
initially-corrected  data  are  input  on  unit  11.  Then  the  following  param¬ 
eters  are  used  in  the  Acceleration,  Filtering,  and  Decimation  Program: 

$SPARMS 

103-12, 

ACC-. FALSE., 

LPF-.TRUE. , 

NBAND-32, 

IDEC-5,  (Other  parameters,  STRTME, 

TMRPLT-.TRUE.,  endtme,  SCLFCT,  PKHT,  TRNC, 

D0PPLT- 1 . ,  etc.,  are  used  as  in  WSPRC.) 

NFFT-64, 

LENWDW-16, 

LAG-4, 

$END 
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Figure  9.  Example  of  Bad  Data  from  24  to  26  Seconds 


This  procedure  will  output  the  lowpass-f iltered  data  on  unit  12.  Next, 
with  the  initially-corrected  data  on  unit  11,  and  the  lowpass-f iltered  data 
on  unit  12,  WSPRC  is  run  with  the  following  parameters. 


$PARMS 

LEVEL=0 , 0 , 1 , 1 , 1 , 
XKNTSP»-.2, 
NAVGP=3 , 

104=0, 

D0PPLT=1 . , 
NFFT=512, 
LENWDW=160, 
LAG=32 , 

$END 


Values  of  SCLFCT,  PKHT ,  and  TRNC  should  be  chosen  to  give  a  readable  plot 
showing  the  spin  lines  (SCLFCT=0. ,  PKHT=4. ,  and  TRNC=1 .  usually  work). 

This  procedure  will  place  the  outputted,  lowpass-f ilter-corrected  data 
on  unit  13.  From  this  stage  in  the  processing, clear  spin  lines  should  be 
visible  on  the  TMR  plot.  An  unsatisfactory  correction  will  result  in  a 
misaligned  central  peak  and  weak  spin  lines.  If  further  correction  is 
needed,  the  above  procedure  should  be  repeated  using  the  recently  corrected 
data  as  input  on  unit  11. 
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3.  INTERPRETING  THE  SPIN  LINES 

Now  Chat  Che  spin  lines  are  clearly  visible  on  the  TMR  plot,  the  next 
step  is  to  identify  them.  The  GSRS  target  has  four  fins,  so  strong  spin 
lines  will  be  seen  at  multiples  of  four  times  the  spin  frequency  fg.  Since 

the  spin  frequency  is  usually  between  8  and  10  Hz,  the  spin  lines  usually  occur 

at  spacings  of  32  to  40  Hz.  Figures  10  through  12  show  examples  of  spin  lines 
at  various  multiples  of  4  fg.  In  all  of  these  TMR  plots  the  spin  lines  occur 
at  intervals  of  about  30  to  40  Hz. 

In  Figure  10  the  -20  fg  spin  line  appears  on  the  positive  side  of  the 
TMR  plot.  This  is  known  as  "foldover."  When  a  spin  line  goes  off  the  TMR  plot 
at  one  side,  it  continues  on  the  ocher  side.  In  Figure  12  spin  lines  appear 
at  multiples  of  up  to  32  fg.  All  of  the  lines  at  multiples  of  20  f  and  higher 
have  been  folded  over.  The  left  and  right  boundaries  of  the  TMR  plot  are 
at  +  PRF/2,  which  in  this  case  is  approximately  +  160  Hz.  So,  if  f  ■  10  Hz, 
the  20  f  spin  line  appears  at  200  Hz.  But  the  plot  goes  to  only  160  Hz, 
so  this  line  will  be  40  Hz  beyond  the  right  edge  of  the  plot,  which  means 
it  will  be  seen  40  Hz  to  the  right  of  the  left  side  of  the  plot,  i.e.,  at 

-120  Hz.  It  is  as  if  the  plot  were  wrapped  on  a  cylinder,  with  the  left  and 

right  sides  joined  together. 

The  spin  lines  are  not  always  straight  and  vertical,  because  the 
spin  frequency  is  not  always  constant.  The  spin  frequency  may  be  in¬ 
creasing  or  decreasing.  If  it  is  increasing,  the  spin  lines  will  be 
spreading  out  and  getting  farther  from  the  central  peak  and  from  each 
other.  If  f  is  decreasing,  the  lines  will  be  getting  closer  to  each 
ocher  and  closer  to  the  center.  This  may  look  confusing  if  there  is 
foldover,  as  in  Figure  12a.  There  the  spin  frequency  increases  from 
the  bottom  to  the  top  of  the  plot.  This  means  that  all  the  positive 
spin  lines  are  slanting  to  the  right.  However,  half  of  the  positive 
spin  lines  are  folded  over,  so  when  they  slant  to  the  right  they  get 
closer  to  the  center,  though  their  frequencies  are  actually  increasing. 

In  Figure  12b  the  spin  frequency  is  decreasing,  so  all  the  positive 
lines  slant  to  the  left,  and  the  negative  lines  slant  to  the  right. 
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Figure  10.  GSRS  Data  from  Site  R-394 


Another  characteristic  to  notice  is  that  the  change  in  spin  frequency 
is  more  pronounced  in  the  higher-order  lines.  For  example,  in  Figure  12a 
the  28  fs  line  slants  strongly  to  the  right,  while  the  4  f  line  looks  al¬ 
most  vertical.  This,  of  course,  is  because  the  lines  are  spaced  at  multiples 
of  4  fs.  Therefore,  when  the  spacing  changes  due  to  a  frequency  change,  the 

28  fg  line  changes  seven  times  as  much  as  the  4  f  line. 
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4.  THE  SPIN  FREQUENCY  SOFTWARE 


Once  Che  spin  lines  have  been  identified,  they  may  be  analyzed  by  the 
WSGSRS  program.  Four  positive  and  four  negative  pairs  of  spin  lines  at  con¬ 
secutive  multiples  of  4  f  are  selected  for  use  in  the  program.  For  example , 

s 

in  Figure  1  the  spin  lines  at  +  4,  +  8,  +  12,  and  +  16  times  the  spin  fre¬ 
quency  should  be  used  by  the  program.  In  Figure  12,  the  +  16,  +  20,  +  24, 
and  +  28  spin  lines  should  be  chosen  because  higher-order  lines  give  more 
accurate  measurements. 

The  WSGSRS  program  tracks  the  peaks  of  the  selected  spin  lines  using 
the  Fourier  transform  software  of  WSPRC  LEVEL(5).  A  noise  level  threshold 
value  THRESH  is  input  in  the  STRCKER  namelist  to  specify  the  power  level  be¬ 
low  which  the  spin  lines  should  be  ignored.  The  proper  value  for  THRESH  is 
determined  by  measuring  Che  height  of  the  noise  on  the  TMR  plot,  and  then 

dividing  this  value  by  the  scale  of  the  plot.  For  example,  in  Figure  10  the 

2 

noise  level  would  be  approximately  THRESH  »  0.04  inches  r  0.03951  in<(v  /Hz) 

2 

-  1.0  v  /Hz.  The  program  determines  for  what  times  the  spin  lines  are  above 

the  specified  threshold,  and  the  spin  frequencies  represented  by  these  spin- 

line  peaks  are  plotted  versus  time.  Figure  13  shows  such  a  spin-frequency 

plot  made  from  the  data  of  Figure  12.  The  jumps  in  frequency  near  19  and 

27.4  seconds  are  false  measurements  caused  by  interference  between  folded 

spin  lines  at  those  times.  It  is  seen  from  the  TMR  plot  in  Figure  12  that 

the  spin  lines  are  superimposed  in  these  regions,  hence  the  spin  frequency 

measurement  will  be  inaccurate  for  these  regions. 

The  software  has  two  methods  of  tracking  the  spin  lines;  automatic  and 

guided.  The  automatic  tracker  tracks  the  spin  lines  at  +  L0WFS  times  the 

spin  frequency  by  searching  for  each  successive  peak  on  the  basis  of  the 

position  of  the  preceding  two  peaks.  The  other  spin-line  peaks  are  searched 

for  at  intervals  of  4  f  from  the  lower-order  lines.  This  automatic  tracker 

s 

requires  distinct  spin  lines  at  +  L0WFS  times  fg.  The  +  4  f  lines  in  Figure  1, 

for  example,  are  suitable  for  the  automatic  tracker. 

Usually  the  spin  lines  at  +  L0WFS  times  f  are  not  completely  clear,  due 

—  s 

to  a  high  noise  level;  or  due  to  interference  by  clutter  or  other  spin  lines. 

In  these  cases  the  automatic  tracker  is  not  used.  Instead,  the  spin-line 
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tracker  is  guided  by  input  data  points  specified  by  the  namelist  arrays 
SPNFRQ  and  TIME,  where  SPNFRQ(I)  equals  f  at  time  TIME(I).  The  program 
uses  straight-line  interpolation  between  these  points  to  estimate  the  spin 
frequency.  The  estimate  is  multiplied  by  L0WFS  to  find  the  starting  point 
for  the  search  of  the  spin  lines.  This  mode  was  used  for  the  data  of 
Figure  12,  and  the  values  used  are  shown  in  Figure  14.  The  values  SPNFRQ 
and  TIME  chosen  were  evaluated  by  measuring  the  frequency  of  the  28  f  line 
in  Figure  12,  and  dividing  by  28.  For  example,  at  26.9  seconds  the  28  f 
line  is  at  -64  Hz.  But  since  the  line  has  been  folded  over,  the  actual 
value  is  -64  +  PRF  =  -64  +  320  *  256,  which  gives  a  spin  frequency  of 
f  ■  256  t  28  ■  9.14  Hz.  Higher-order  spin  lines  are  used  for  such  cal¬ 
culations  because  they  give  more  accurate  measurements.  Enough  of  these 
input  data  points  should  be  chosen  so  that  the  actual  spin  frequency  stays 
within  about  0.1  Hz  of  the  straight-line  interpolated  values.  One  must  have 
T1ME(1)  <  STRTME  and  TIME (last)  £  ENDTME,  since  the  program  does  not  extra¬ 
polate. 

Following  is  a  table  describing  the  use  of  the  input  parameters  by  the 
spin-line  tracking  software. 

NAMELIST /TRCKER/  Parameters: 

Description 

Determines  whether  automatic  tracking  is  used. 

If  .TRUE.,  the  program  tracks  the  spin  lines 
by  following  the  peaks  of  the  lowest  spin 
line  (L0WFS) .  But  this  automatic  tracker 
is  unpredictable  if  the  data  are  bad  or 
if  the  spin  lines  are  ambiguous.  For 
such  cases  use  AUT0T=. FALSE.  When 
AUT0T=. FALSE. ,  breakpoints  specifying 
time  and  spin  frequency  are  read  in  through 
the  variables  TIME  and  SPNFRQ.  Straight- 
line  interpolation  between  these  points 
then  guides  the  tracker. 


Name  (Type,  Default) 
$TRCKER 

AUT0T  (LOG. , .TRUE.) 
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Ln  ^  w  to  H 


Actual  Spin  Frequency 


Input  Spin  Data  Points  H - f- 

and  Interpolated  Lines 


Time  (sec) 


TIME (I) 

SPNFRQ( 

18.2 

8.68 

20.8 

9.39 

23.0 

9.75 

24.6 

9.85 

25.05 

9.79 

26.9 

9.14 

Figure  14.  Actual  and  Estimated  Spin  Frequency 
for  Data  of  Figure  12 


Name  (Type,  Default) 


Descript  ion 


L0WFS(INT. ,4) 


The  multiple  of  the  spin  frequency  represented 
by  the  lowest  frequency  line  to  be  tracked. 

The  tracker  follows  lines  at  i(L0WFS+O, 

+4, +8, +12);  8  lines  in  all. 


STRTHZ (REAL , 8 . ) 


1098 (INT. .98) 
THRES (REAL, 100.) 


TIME (REAL, 20*-l.) 


The  spin  frequency  (Hz)  at  the  start  of  the  trac 
ing  period.  The  auto-tracker  starts  looking  for 

spin  lines  at  tSTRTHZ*L0WFS .  If  AUT0T= 

.FALSE.,  STRTHZ  need  not  be  specified,  as  it 
is  calculated  from  the  input  breakpoints. 

The  temporary  unit  to  which  the  tracking  data 
are  written. 

2 

The  noise  level  threshold  (v  / Hz)  of  the  Doppler 
peaks  power  spectrum.  When  the  spin-line  power 
falls  below  this  level  it  is  ignored.  The  value 
for  this  parameter  must  be  carefully  chosen  from 
the  previous  TMR  plots. 

The  array  containing  the  breakpoint  times  to 
be  used  in  calculating  the  spin  frequency  for 
the  tracker.  The  spin  frequency  is  SPNFRQ(I) 
at  time  TIME(I).  Only  the  values  of 
TIME  preceding  the  first  negative  value  are 
used. 


SPNFRQ (REAL, 20*8. )  The  spin  frequencies  (Hz)  at  the  times 

TIME.  These  valuer  must  be  estimated 
from  the  previous  TMR  plots. 

PLTLNS (LOG ., .FALSE . )  If  PLTLNS= .TRUE . ,  a  frequency  plot  and 

a  power-level  plot  will  be  generated 
for  each  of  the  eight  spin  lines. 


Since  the  WSGSRS  software  uses  WSPRC  LEVEL (5)  software,  the  LEVEL (5) 
namelist  $PARMS  is  required  as  input.  The  deck  set  up  for  input  to  WSGSRS 
is  as  follows; 
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@  E^F  Card 
Header  Card 
$PARMS 

LEVEL-0,0,0,0,1, 

STRTME=. . . ,ENDTME=. . . ,103=. . .  , 

106=0, 

TMRPLT=  .FALSE.,  (We  assume  that  suitable  TXR  plots 
D0PPLT=1.  ,  have  been  previously  generated.) 

NFFT=. . .  ,  LENVDW=. . . ,  LAG=.  .  .  , 

SEND 

STRCKER 

Tracker  parameters 

SEND 

@E0F  or  other  master-space  card 

(Usual  values  for  the  LEVEL(5)  parameters  are:  NFFT=512,  LENVDW=160,  LAG=32.) 

Figure  15  shows  the  printer  output  from  the  analysis  of  the  data  of 
Figure  12.  The  WSPRC  parameters  are  printed,  then  the  WSGSRS  parameters, 
then  the  data  from  the  spectral  tracking.  These  tracking  data  are  written 
to  the  unit  specified  by  the  variable  1098,  and  may  be  plotted  by  specifying 
PLTLNS=.TRUE.  The  data  are  written  in  groups  of  17  values  in  the  order 

time, 

freq_L’  ^P-l*  freqL’  “PL*  “P-O+A) . 

freq(L+8)*  “^+8)  ’  f req-(L+12)  ’  amp-(L+12)  *  freq(L+12)’  “Pa+lZ)’ 

where  the  subscript  denotes  the  spin-frequency  multiple,  and  L=L0WFS. 

Next,  the  intervals  for  which  the  spin  line  amplitudes  exceed  the  noise  level 
are  listed.  For  each  of  the  8  spin  lines,  20  pairs  of  times  are  given.  A 
value  of  -1  means  that  no  more  intervals  were  found. 

A  spin -frequency  tape  is  also  output  by  the  WSGSRS  program.  The  tape, 
on  unit  15,  contains  an  80-character  header,  followed  by  time/spin-frequency 
data  pairs  in  the  (1X,2E20.12)  format.  The  spin-frequency  data  written  out 
are  median  averages  which  are  superimposed  on  the  spin-frequency  plot.  A 
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Figure  13a.  Output  from  GSRS  Analysis  of  Data  of  Figure  12 
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Figure  15d .  Output  from  GSKS  Analysis  of  Data  of  Figure  12 
(Continued) 


value  of  spin  frequency  is  written  for  each  data  time  plotted.  If  an  even 
number  of  spin-frequency  values  is  plotted  at  a  particular  time,  then  the 
average  of  the  interior  two  values  is  written  to  tape.  If  an  odd  number  of 
spin-frequency  values  is  plotted,  then  the  center  (median)  value  is  written 
to  tape. 
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5.  THE  CONING-MOTION  SOFTWARE 


The  TCM  software  makes  target  coning  measurements  via  the  range  residual 
method.  The  range  residual  is  the  phase-derived  range  deviation  from  zero 
left  over  after  the  c.g.  trajectory  motion  is  subtracted  from  the  phase- 
derived  range  data.  The  instantaneous  range  residual  AR  gives  the  instan¬ 
taneous  coning  angle  according  to  the  relationship 

AR  -  l  sinft  sinft 
c  o 

where 

AR  *  range  residual, 

l  •  length  from  c.g.  to  base  of  missile, 

*  aspect  angle,  and 

*  the  projection  of  the  target  coning  angle  into  the  plane  formed 
by  the  radar  line  of  sight  and  the  trajectory  velocity  vector. 


In  this  software  the  trajectory  of  the  c.g.  is  determined  by  a  polynomial  fit 
to  the  square  of  the  range.  The  TCM  software  is  divided  into  two  main  pro¬ 
grams.  The  first  is  CRSQRD,  which  performs  the  standard  range-residual 
measurement.  The  second  is  SLR,  which  performs  several  range-squared  fits 
in  a  process  designed  to  minimize  fitting  error. 


Three  data  files  are  used  as  input  by  the  TCM  program:  (1)  the  WSMR  site 
radar  data  tape;  (2)  the  lowpass  filtered  raw  amp-phase  data  file;  and  (3)  the 
aspect-angle  data  file.  The  use  of  the  aspect-angle  file  is  optional  when  the 
CRSQRD  program  is  run. 


5.1.1  WSMR  Site  Tape 

The  WSMR  site  radar  data  tape  is  assigned  to  unit  1.  The  program  uses 
the  range  data  from  this  tape  to  make  an  initial  trajectory  fit. 


5.1.2  LPFRAW  File 

The  lowpass-f iltered  raw  amp-phase  (LPFRAW)  file  is  assigned  to  unit  101 
(normally  unit  11).  The  data  from  this  file  are  used  to  refine  the  trajectory 
fit.  The  fitting  procedure  requires  amp-phase  data  which  represent  the  un¬ 
corrected  (raw)  base  return  only,  without  interference  from  other  scattering 


45 


centers,  such  as  the  spinning  fins.  The  creation  of  this  desired  data  file 
is  accomplished  in  three  steps:  (1)  perform  the  initial  correction  of  the 
data  using  WSPRC;  (2)  lowpass  filter  the  corrected  data  so  that  all  but  the 
base  return  is  filtered  out;  and  (3)  use  WSPRC  to  reverse  the  trajectory 
correction  done  in  step  (1)  so  that  the  data  are  again  "raw"  (i.e.,  so  that 
all  trajectory  compensation  is  undone). 

The  first  of  the  three  steps  is  done  exactly  as  previously  described  in 
Section  2.1.  This  step  is  called  the  initial  correction.  The  cumulative 
spline  breakpoints  are  output  by  WSPRC  on  unit  16,  and  must  be  saved  for  use 
in  step  (3).  If  two  or  more  splines  are  necessary  for  the  initial  correction, 
they  all  must  be  saved.  The  initial  correction  is  satisfactory  when  the  cen¬ 
tral  peak  cn  the  TMR  plot  stays  between  -20  and  +20  Hz. 

In  the  second  step  this  initial -correction  file  is  lowpass  filtered  by 
using  it  as  input  to  the  Acceleration,  Filtering,  and  Decimation  Program 
with  the  following  parameters: 

$PARMS 

ACC-. FALSE. , 

LPF-.TRUE. , 

TMRPLT- . TRUE . , 

NBAND-32, 

ID EC-1, 

NFFT-256 , 

LENWDW-64, 

LAG-16, 

HRZ1— 80.  .HRZ2-80. , 

$END 

The  third  step  is  to  use  the  output  lowpass-filtered  file  from  the 
second  step  as  input  to  LEVEL (4)  of  WSPRC  in  order  to  apply  the  reverse 
of  the  splines  saved  from  step  (1).  If  two  splines  were  used,  they  are 
assigned  to  units  1(34  and  1(35.  If  only  one  spline  was  used,  it  is  as¬ 
signed  to  unit  1(35,  and  parameter  104-0  is  set.  The  input  amp-phase  file 
is  assigned  to  unit  101.  The  following  parameters  would  be  used  to  reverse 
the  correction  applied  by  the  spline  on  unit  105: 
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LEVEL  -0,0, 0,1,0, 
104-0, 

C0EFS-O. ,-l. ,-l. , 
SEND 


The  values  of  C0EFS  indicate  chat  the  negative  of  the  spline  is  to  be  used. 

The  outputted  LPFRAW  file  is  saved  on  103.  This  is  the  file  to  be  used 
as  input  to  the  TCM  program. 

5.1.3  Aspect-Angle  Data  File 

The  aspect-angle  data  file  is  assigned  to  unit  104  (normally  unit  14) 
The  aspect  angles  are  used  to  convert  the  range  residual  into  coning  angle 
according  to  the  relation 

n  ■  Sin  ^[AR/Usiafl  )]  . 
c  o 

The  aspect-angle  file  consists  of  a  96-character  header  followed  by 
time— aspect  angle  pairs  written  in  the  following  manner: 

DIMENSI0N  HEADER(16) 

REWIND  104 
WRITE(I04)  HEADER 
1  T  -  new  time 

A  -  new  aspect  angle 
WRITE (104)  T.A 
IF  (more  data)  G0  T0  1 
CALL  ENDFYL(I04) 

5.2  THE  CRSQRD  PROGRAM 

Parameters  are  input  to  the  program  CRSQRD  through  namelist  $PARMS.  A 
header  card  precedes  the  namelist.  The  namelist  parameters  and  an  explanation 
of  their  use  are  as  follows. 


Name  (Mode,  Default) 


Description 


STRTME (REAL , - 1 . E30 

ENDTME(R£AL, 1.E30) 

KJ1CINT.11) 

KJ3CINT.13) 

NAVGR(INT,40) 

NAVGP(INT.O) 

I04(INT,14) 

CGL(REAL,1.5) 

PLOT  4 (L0G , .TRUE. ) 

XI (REAL, 0.) 

X2CREAL.0.) 


The  biased  first  time  (seconds)  of  the  time 
interval  to  be  processed.  Defaults  to  the 
first  data  time  available.  The  time  bias 
is  obtained  from  Che  LPFRAW  data  described 
in  Section  2.2. 


The  biased  last  time  (seconds)  of  the  time  in¬ 
terval  to  be  processed.  Defaults  to  the  last 
data  time  available. 

The  input  unit  assigned  to  the  LPFRAW 
data. 

The  output  unit  assigned  to  the  output 
corrected  amp-phase  data. 

The  number  of  range  data  points  averaged 
into  each  group  for  the  initial  range- 
squared  fit. 

The  number  of  phase  data  points  averaged 
into  each  group  for  the  refined  range- 
squared  fit.  If  not  positive,  defaults 
to  NAVGR/5 . 

The  input  unit  assigned  to  the  aspect- 
angle  data.  If  104*0 ,  the  plot  of 
angular  variation  is  not  made. 

The  length  in  meters  from  the  c.g.  to  the 
base  of  the  missile. 

If  PLOT 4*. FALSE. ,  the  first  three  output 
plots  will  not  be  made. 

If  X1#0.,  the  left  boundary  of  the  output 
plots  will  be  XI  or  the  first  output  data 
time,  whichever  is  smaller.  If  X1*0., 
the  boundary  is  determined  by  the  program. 

If  X2^J.,  the  right  boundary  of  the  out¬ 
put  plots  will  be  X2  or  the  last  output 
data  time,  whichever  is  greater.  If  X2*0., 
the  boundary  is  determined  by  the  program. 
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The  output  from  CRSQRD  consists  of  printout,  four  plots,  and  a  file  on 
unit  103.  This  output  file  is  a  scratch  file  containing  amp-phase  data  cor¬ 
rected  by  the  range-squared  fit  computed  by  the  program. 

The  printout  lists  the  parameters  of  the  polynomials  of  the  initial  and 
updated  range-squared  fits.  The  input  parameters  and  the  identities  of  the 
input  and  output  files  are  also  listed. 

Figures  16a,  b,  c,  and  d  are  examples  of  the  plots  which  program  CRSQRD 
produces  to  show  the  steps  used  in  the  measurement  of  coning  motion.  The 
first  plot  (Figure  16a)  shows  two  curves  representing  the  range  residual 
due  to  the  initial  range-squared  trajectory  fit.  One  curve  is  the  residual 
in  the  range  data;  the  other  is  the  residual  in  the  phase  data.  The  phase 
residual  is  the  smoother  curve. 

The  phase  residual  is  added  to  the  initial  fit  to  provide  refined  data 
for  another  range-squared  polynomial  fit.  The  residual  from  this  updated 
fit  is  plotted  in  the  second  plot  (Figure  16b).  The  third  plot  (Figure  16c) 
shows  the  residual  range  variation  of  the  base  when  the  new  fit  is  applied 
to  the  phase  data  without  grouping  and  averaging  of  the  data. 

If  aspect-angle  data  are  available  to  the  program  (i.e. ,  if  I04#»  ,  then 
the  range  variations  of  the  base  as  shown  in  the  third  plot  are  converted  to 
angular  motion,  and  are  plotted  in  the  fourth  plot  (Figure  16d).  This  plot 
represents  the  desired  output  from  the  program.  The  first  three  plots  are 
useful  primarily  for  ensuring  that  the  fitting  procedure  was  successful. 

5.3  CHARACTERISTICS  OF  RANGE-SQUARED  FITTING 

In  order  for  targec  coning  motion  to  be  measured,  it  is  crucial  that 
the  c.g.  trajectory  be  accurately  determined.  The  generalized  trajectory 
equation  for  a  constant  acceleration  vector  (e.g.,  for  exo-atmospheric 
free  fall  over  a  flat  earth)  is 

R  =  Jk  +  Bt  +  Ct2  +  Dt3  +  Et4  ,  (1) 

where  R  is  the  range  from  a  reference  point  (e.g.,  from  a  radar  site).  With¬ 
in  the  atmosphere,  however,  this  equation  is  good  only  as  an  approximation 
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to  the  trajectory  over  short  time  intervals.  This  is  because  of  the  non¬ 
constant  acceleration  due  to  atmospheric  friction.  Over  a  short  time  in¬ 
terval,  though,  one  may  make  the  assumption  that  the  acceleration  is  nearly 
constant,  and  hence,  the  trajectory  may  be  approximated  by  fitting  a  fourth- 
order  polynomial  to  the  square  of  the  range. 

When  the  root-polynomial  fit  is  subtracted  from  the  range  data,  the  range 
residual  remains.  If  the  trajectory  fit  were  perfect,  and  the  range  measure¬ 
ments  'jere  to  the  base  of  the  missile,  this  residual  would  correspond  exactly 
to  the  target  coning  motion  according  to  the  relationship 


ft  =  sin  "H  AIV  (2.  sinft  )] 
c  eg  o 

where 

AR  =  range  residual  » 

£,  =  distance  from  target  c.g.  to  base, 

eg 

ft  =  aspect  angle, 
o 

ft  =  projected  coning  angle. 


If  the  time  interval  used  for  fitting  were  too  short,  the  polynomial  would 
follow  the  coning  motion  in  addition  to  the  trajectory  motion,  so  the  resi¬ 
dual  would  be  less  than  the  actual  coning  motion.  Conversely,  if  the  time 
interval  were  too  long,  the  polynomial  would  not  be  able  to  completely  com¬ 
pensate  for  the  trajectory  (due  to  a  variable  acceleration)  so  the  residual 
would  be  greater  than  the  actual  coning  motion.  The  problem,  then,  is  to 
decide  when  the  polynomial  fit  correctly  compensates  for  the  trajectory. 

The  polynomial  being  used  for  the  fit  is  of  fourth  order,  so  it  will  have 
five  extrema  as  illustrated  in  Figure  17.  The  polynomial  can  be  represented 


by  a  linear  combination  of  the  zeroth  through  fourth-order  Legendre  poly¬ 
nomials  (Pq  through  P^) .  The  Legendre  polynomials  are  orthogonal  polynomials 
defined  within  a  finite  interval.  After  these  five  Legendre  polynomials  are 


removed  from  the  trajectory,  the  most  significant  term  remaining  will  be  the 

3  5 

fifth-order  Legendre  polynomial,  P,.  ■  1/8  [15x  -  70x  +  63x  ],  shown  in 

Figure  18.  The  range  residual  should  have  this  six-extrema  shape  when  the 
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Figure  17.  Fourth-order  Polynomial 


Figure  18.  Fifth-Order  Legendre  Polynomial 


fitting  interval  is  sufficiently  long  so  that  the  fitting  polynomial  does 
not  remove  all  of  the  trajectory  motion  (but  when  the  fitting  interval  is 
sufficiently  short  for  the  term  to  still  be  negligible).  Figure  19  is 
an  example  of  such  a  residual.  In  that  example  the  time  interval  for  fit¬ 
ting  was  too  long,  and  the  residual  has  the  characteristic  shape  of 
Figure  18..  The  time  interval  could  be  shortened  to  improve  the  trajectory 
fit,  but  there  is  a  limit  to  how  short  the  interval  can  be.  This  limit  is 
reached  when  the  fitting  polynomial  is  able  to  follow  the  target  coning  mo¬ 
tion  in  addition  to  fitting  the  trajectory  motion.  The  polynomial  in 
Figure  17  has  two  cycles  of  an  oscillation.  Therefore,  if  the  time  inter¬ 
val  corresponds  to  two  precession  cycles,  the  precession  motion  will  be  in¬ 
cluded  in  the  trajectory  compensation  and  will  not  appear  in  the  residual. 
In  fact,  since  the  fifth-order  residual  has  two  and  one-half  cycles,  a 
fitting  interval  must  be  longer  than  two  and  one-half  cycles  of  a  target 
oscillation  in  order  to  resolve  that  oscillation.  For  example,  if  the  pre¬ 
cession  period  were  3  seconds,  a  fitting  interval  of  at  least  7.3  seconds 
would  be  required  to  separate  the  coning  motion  from  the  trajectory  motion 
in  the  measured  range. 

The  nutation  of  the  GSRS  missile  is  so  fast  that  there  is  no  problem 
in  separating  it  from  the  trajectory  motion.  The  nutation  period  is 
usually  about  0.4  sec,  so  a  fitting  interval  of  several  seconds  is  adequate 
to  separate  it  from  the  trajectory  motion,  because  then  many  nutation  cy¬ 
cles  are  contained  within  one  fitting  interval.  Figure  20  shows  nutation 
motion  over  a  six-second  interval. 

A  longer  fitting  interval  is  necessary  to  resolve  precession  motion. 
Figure  16d  shows  a  13-second  interval  which  resolves  precession  motion  with 
a  period  of  about  3  seconds.  However,  longer  precession  periods  can  make 
it  impossible  to  resolve  the  precession  motion  so  simply.  This  happens 
when  the  trajectory  deviates  from  the  model  of  Eq.  (1)  to  the  extent  that 
the  polynomial  cannot  correctly  fit  the  trajectory  over  a  time  interval 
which  is  longer  than  2  1/2  precession  cycles.  This  problem  is  illustrated 
in  Figures21a,  b,  and  c.  Shorter  fitting  intervals  are  used  in  each  suc¬ 
cessive  plot  to  remove  the  characteristic  six-extrema  residual  shape.  This 
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Figure  19.  Residual  from  a  10-second  Fitting  Interval 
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Figure  20.  Nutation  Motion 
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Figure  21a.  Residual  from  a  7-second  Fitting  Interval 
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Figure  21b.  Residual  from  a  6-second  Fitting  Interval 
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Figure  21c.  Residual  from  a  5-second  Fitting  Interval 
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residual  shape  finally  disappears  in  Figure  21c,  for  which  Che  fitting  in¬ 
terval  has  been  reduced  to  5  seconds.  But  in  Figure  21c  the  precession 
motion  is  obviously  removed  also.  In  Figures  21a  and  21b  2-1/2  cycles  of 
motion  are  always  seen,  so  it  is  impossible  to  determine  which  cycles  are 
due  to  precession  and  which  are  due  to  trajectory  residual. 

5.4  THE  SLR  PROGRAM 

There  is  a  solution  to  the  problem  of  separating  trajectory  and 
coning  motion  components.  As  explained  earlier,  when  the  fitting  inter¬ 
val  is  too  long  and  the  polynomial  cannot  follow  all  the  trajectory  motion, 
the  main  feature  of  the  residual  will  be  a  shape  similar  to  that  of  the 
Legendre  polynomial  P^.  It  is  also  known  that  the  polynomial  has  five 
roots.  At  each  of  these  roots  the  error  in  the  range  residual  represented 
by  P^  will  be  small.  In  other  words,  the  systematic  fitting  error  will  be 
a  minimum  at  the  locations  of  these  roots.  Therefore,  the  values  of  the 
range  residuals  at  the  locations  of  these  roots  should  be  good  estimates 
of  the  true  range  residual  due  to  target  coning.  At  the  first  and  last 
roots,  the  slope  of  the  P^  polynomial  is  very  large,  so  that  a  slight  er¬ 
ror  in  locating  the  root  could  give  a  large  error  in  the  corresponding 
range  residual.  For  this  reason,  the  range  residuals  at  these  roots  are 
not  used  as  estimates.  The  residuals  at  the  three  interior  roots  give 
more  trustworthy  estimates. 

The  three  roots  of  the  P^  polynomial  are  easily  calculated.  They  lie 
at  fractions  0.23077,  0.5,  and  0.76924  of  the  total  fitting  interval. 

Figure  22  shows  the  locations  of  the  roots  marked  on  the  data  from  Figure  19. 
If  this  fitting  interval  is  shifted  to  the  left  or  right,  a  new  fit  can  be 
made,  yielding  three  more  range  residual  estimates  at  new  times.  In  this  man¬ 
ner,  the  entire  precession  motion  can  be  constructed  by  sliding  a  fitting  in¬ 
terval  along,  and  plotting  the  three  range  residual  estimates  of  each  fit 
along  the  way.  Figure  23  is  just  such  a  plot.  Six  fits  were  made  to  the 
data  at  starting-time  increments  of  0.42  seconds.  Each  fitting  interval  was 
10  seconds  in  length.  Three  small  sections  of  data  were  plotted  from  each 
fitting  interval.  These  sections  of  data  were  centered  on  the  Legendre  roots 
of  the  various  fits.  A  line  has  also  been  plotted  which  passes  through  the 
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Figure  22.  Roots  of  Legendre  Polynomial  P5  Marked 
on  Residual  Data 
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average  value  of  each  data  section.  This  line  represents  the  desired  esti¬ 
mate  of  the  precession  motion.  This  method  of  target  coning— mot  ion  measure¬ 
ment  is  called  the  Sliding  Legendre  Roots  (SLR)  method. 


The  plots  of  Figures  21  and  23  were  made  from  the  same  data.  Compar¬ 
ison  of  these  plots  illustrates  the  value  of  the  SLR  method  in  resolving 
precession  motion. 

The  SLR  program  is  executed  by  the  TCM  software  when  parameter 
S0LAR=.TRUE.  is  set  in  namelist  $PARMS.  Programs  CRSQRD  and  SLR  have  the 
following  parameters  in  common:  101,  104,  NAVGR,  NAVGP ,  CGL,  XI, and  X2. 
Program  SLR  does  not  use  parameters  STRTME  and  ENDTME.  Instead,  it  uses 
several  other  parameters,  some  of  which  are  illustrated  in  Figure  24.  A 
complete  list  of  the  parameters  used  by  program  SLR  follows. 


Name  (Mode,  Default)  Description 

S0LAR(L0G, .FALSE.)  If  S0LAR= . TRUE . ,  program  SLR  is 

called  to  perform  Sliding  Legendre 
Roots  processing.  If  S0LAR=. FALSE. , 
program  CRSQRD  is  called  to  perform 
standard  range-squared  fit  analysis. 


IF  S0LAR= . TRUE . ,  the  following  parameters  are  used: 


1 01 (INT, 11) 
I04(INT,14) 
NAVGR(INT,40) 
NAVGP ( INT, 0) 
CGL (REAL, 1.5) 
XI  (REAL,0.) 
X2  (REAL,0.) 


Same  as  in  CRSQRD  (see  Section  5.2) 


TM1(REAL,0. ) 


TM2 (REAL,0. ) 


Value  of  STRTME  for  first  range- 
squared  iteration.  TM1  must  be 
positive. 

Value  of  STRTME  for  last  range- 
squared  iteration.  Defaults  to 
TM1  +  .25*TMLEN*(1.  -  l./ITR). 


ITR(INT,6) 


Number  of  range-squared  iterations 
to  be  performed.  Maximum  is  20. 


TMLEN (REAL , 10 . ) 


The  length  in  seconds  of  each  range- 
squared  fit. 
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TB0UND (REAL, 0.15) 


Data  within  +  TB0UND  seconds  of  each 
root  is  plotted.  Defaults  to  0.15  or 
(TM2-TM1)/ (2.*(ITR-1)) ,  whichever  is 
smaller. 

SPLT (LOG,. TRUE.)  If  SPLT=.TRUE.,  the  sections  of  data 

centered  at  the  roots  are  plotted. 

LPLT (LOG. .TRUE. )  If  LPLT=.TRUE.,  line  segments  are 

plotted  joining  the  averaged  centers 
of  the  sections  of  data. 

A  desired  SLR  plot,  like  the  one  in  Figure  23,  can  usually  be  made  by 
setting  S0LAR= .TRUE. ,  defining  TM1,  and  letting  all  the  other  parameters 
default.  If  data  are  not  available  up  to  TM2  +  TMLEN  seconds,  processing 
will  stop  short  of  completion.  It  is  important  that  a  suitable  length  of 
time  be  chosen  for  the  fitting  interval  TMLEN.  In  Figure  23  the  default 
value  of  10  seconds  was  used.  Due  to  considerations  previously  discussed, 
TMLEN  should  be  long  enough  so  that  the  P,.  polynomial  is  the  dominant  shape 
in  the  residual.  Usually  for  the  GSRS,  10  seconds  or  longer  will  be  suitable, 
and  the  fitting  interval  will  contain  approximately  4  or  more  precession 
periods. 

5.5  EFFECTIVE  USE  OF  THE  CONING-MOTION  SOFTWARE 

The  CRSQRD  and  SLR  programs  are  an  eff§ctive  combination  for  measuring 
target  coning  motion.  Program  CRSQRD  should  be  used  whenever  possible  to  ac¬ 
curately  measure  coning  motion.  The  fitting  interval  used  should  be  the 
longest  interval  possible  which  does  not  introduce  a  systematic  error  into 
the  trajectory  compensation.  The  nutation  of  the  GSRS  missile  is  invariably 
rapid  enough  to  be  easily  resolved  by  the  CRSQRD  method.  In  some  cases, 
however,  it  will  be  impossible  for  the  CRSQRD  program  to  resolve  the  slower 
precession  motion.  In  these  cases,  program  SLR  should  be  used  to  resolve 
the  slower  motion. 

Whenever  possible,  concurrent  data  from  several  radars  should  be  ana¬ 
lyzed  and  compared  to  confirm  observations  and  to  improve  accuracy.  It 
should  be  remembered  that  each  radar  sees  the  projection  of  the  coning  motion 
into  a  different  plane,  so  the  results  will  not  correspond  identically.  The 


general  pattern  of  the  coning  motion,  however,  should  be  the  same  in  the 
analysis  of  each  radar's  data. 
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6.  SOFTWARE  ON  THE  RADAR  GRAPHICS  LABORATORY 


Programs  WSGSRS  and  TCM  have  been  implemented  as  tasks  on  the  PDP-11/55 
Radar  Graphics  Laboratory  (RGL) .  The  usage  of  these  tasks  on  the  RGL  and 
on  the  Univac  1108  is  identical  except  for  a  few  features  which  will  now  be 
described. 

The  tasks  should  be  run  from  the  graphics  terminal  (as  opposed  to  the 
printer  terminal) .  The  software  accepts  the  input  parameters  from  the  ter¬ 
minal  by  displaying  the  parameter  name  followed  by  its  present  value.  A 
carriage  return  on  the  keyboard  leaves  the  parameter  unchanged.  To  expli¬ 
citly  set  the  parameter,  type  the  desired  value  followed  by  a  carriage  re¬ 
turn.  After  all  parameters  in  the  namelist  have  been  entered  in  this  manner, 
the  procedure  may  be  repeated  to  correct  any  mistakes  which  may  have  been 
made  in  typing  in  the  values.  Headers  and  filenames  are  also  entered  in 
this  manner.  On  the  RGL  files  are  specified  by  file  names  during  execution 
of  the  software,  rather  than  by  logical  unit  assignments  specified  in  control 
cards. 

The  parameters  PAUSE  and  HRDCPY  have  been  added  to  the  namelists  of 
these  tasks  on  the  RGL.  The  logical  variable  PAUSE  determines  whether  or 
not  the  software  will  pause  at  the  end  of  each  plot.  If  PAUSE  equals  .TRUE., 
execution  of  the  task  will  pause  at  the  completion  of  each  plot  page,  pending 
the  receipt  of  any  keystroke  from  the  terminal,  upon  which  execution  resumes. 
Logical  variable  HRDCPY  determines  whether  or  not  automatic  hard  copying  of 
plots  is  invoked.  The  default  values  are  PAUSE*. FALSE,  and  HRDCPY* . TRUE . , 
so  chat  in  the  default  case  plot  pages  are  automatically  copied  and  no  action 
is  required  of  the  operator. 

6.1  THE  SPIN- FREQUENCY  TASK 

The  spin -frequency  task  WSGSRS  is  executed  on  the  RGL  by  the  following 
command : 

@RG : [ 100 , 100 ] RUN  SPIN 

This  software  uses  the  following  input  parameters  which  are  not  implemented 
in  the  Univac  version. 
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Name  (Mode,  Default) 
YMIN(REAL,6.) 


Description 

The  lower  Limit  (in  Hz)  of  the  spin- 
frequency  plot. 

YMAX(REAL, 12 . )  The  upper  limit  (in  Hz)  of  the  spin- 

frequency  plot. 

PRNTSP (L0G. , .FALSE. )  Spin-frequency  data  from  the  scratch 

file  are  not  printed  unless  PRNTSP 
equals  .TRUE. 

The  LEVEL (5)  namelist  $PARMS  on  the  RGL  contains  only  parameters  rele¬ 
vant  to  spin- frequency  processing.  These  are  STRTME,  ENDTME,  NFFT,  LENWDW, 
LAG,  C0HRNT,  ZDSUP,  IWGHT,  HR21,  HRZ2 ,  XLEN,  YLEN,  PAUSE,  and  HRDCPT . 

The  spin-frequency  software  asks  whether  or  not  the  averaged-spin  file 
is  to  be  created.  If  the  answer  is  yes,  the  output  file  is  created  with  a 
name  given  by  the  operator. 

6.2  THE  CONING-MOTION  TASK 

The  coning-motion  task  TCM  is  executed  on  the  RGL  by  the  following 
command > 

@RG: [100, 100] RUN  TCM 

One  additional  input  parameter  is  used  in  this  version  of  program  TCM.  The 
logical  parameter  CTAPE  allows  for  the  creation  of  an  unformatted  data  file 
consisting  of  a  100-character  header  followed  by  time/coning  angle  data  pairs. 
The  data  are  the  points  which  are  plotted  in  the  angular  variation  plot.  If 
CTAPE  equals  .TRUE.,  the  data  file  is  created  on  the  system  disk  with  the 
name  F0RO1O.DAT. 

6.3  THE  SPLINE  FILE  TASK 

Task  SPFILE  has  been  implemented  on  the  RGL  to  perform  the  creation  and 
editing  of  spline  files  by  hand.  This  task  is  executed  by  the  following 
command  : 

@RG : [100, 100] RUN  SPLN 

The  task  asks  for  spline  data  which  it  writes  to  a  spline  file  in  the  WSPRC 
format.  The  default  data  values  presented  by  the  program  come  from  an 


existing  spline  file  (in  the  case  of  file  editing),  or  from  calculations 
(in  the  case  of  file  creation).  When  editing  a  spline  file,  breakpoints 
may  be  inserted,  deleted,  and  altered.  The  program  asks  for  the  breakpoint 
time,  and  then  for  the  R  value.  If  a  new  file  is  being  created,  the  pro¬ 
gram  displays  the  current  default  range  value  calculated  according  to 
the  relation 

RR'Ri  +  i(\  +  V«R-V 

where  the  subscript  R  denotes  the  new  breakpoint,  and  the  subscript  L  de¬ 
notes  the  last  breakpoint.  This  feature  simplifies  the  creation  of  hand 
splines  because  all  that  is  required  is  to  read  time  and  R  from  the  Doppler 
peaks  plot,  to  enter  these  values,  and  to  use  the  default  values  for  range. 
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7.  MEASUREMENT  ACCURACIES 


Spin-frequency  and  coning-mot ion  measurements  of  GSRS  motion  have  not 
been  made  during  the  boost  phase  of  the  missile's  flight.  This  is  because 
during  the  boost  phase  the  radar  skin  return  is  obscured  by  the  exhaust 
gases.  Accurate  measurements  can  be  made,  however,  after  the  boost  phase. 

7.1  SPIN  FREQUENCY 

Figure  13  shows  the  spin  frequency  measurement  of  a  GSRS  flight.  The 
accuracy  obtainable  is  on  the  order  of  +0.03  Hz,  which  is  the  average  width 
of  the  composite  line  formed  by  superimposing  the  several  individual  spin- 
frequency  measurements  made  by  the  WSGSRS  program. 

7.2  CONING  MOTION 

Figures  16d  and  20  show  good  measurements  of  nutation  motion.  Some 
modulation  will  be  seen  in  the  amplitude  of  this  motion  since  the  angle 
measured  is  the  projection  of  the  coning  angle  into  the  plane  formed  by 
the  radar  line  of  sight  and  the  GSRS  trajectory  vector.  The  accuracy  of 
the  measurement  of  this  projected  angle  is  on  the  order  of  +  0.3°  (-5  mil) 
which  is  the  approximate  amplitude  of  the  noise  seen  in  the  two  plots.  The 
measurement  of  the  slower  precession  motion  is  less  accurate.  The  estimated 
accuracy  of  precession  measurement  in  Figure  23  is  +  0.5°. 
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